Nonlinear microwave response of MgB 2 
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We calculate the intrinsic nonlinear microwave response of the two gap superconductor MgB2 in 
the clean and dirty limits. Due to the small value of the tv band gap, the nonlinear response at low 
temperatures is larger than for a single gap Bardeen-Cooper-Schrieffer (BCS) s-wave superconductor 
with a transition temperature of 40 K. Comparing this result with the intrinsic nonlinear d-wave 
response of YBa2Cu3C>7 (YBCO) we find a comparable response at temperatures around 20 K. Due 
to its two gap nature, impurity scattering in MgB2 can be used to reduce the nonlinear response if 
the scattering rate in the it band is made larger than the one in the a band. 

PACS numbers: 74.25.Nf, 74.70.Ad, 84.40.Dc 



High-T c cuprate thin films are being used for construct- 
ing microstrip resonators and filters in the microwave 
regime. One of the limiting factors is their nonlinear 
response leading to undesirable harmonic generation and 
intermodulation. While this nonlinear response is often 
associated with weak links, there is an intrinsic nonlincar- 
ity which sets a lower limit on what can be achieved. This 
has previously been discussed for the cuprates where the 
d-wave nature of the gap leads to an increase in the intrin- 
sic nonlinear response at low temperatures 0,0- The re- 
cent discovery of s-wave superconductivity in MgB2 with 
a comparatively high critical temperature of T c = 40 K 
and progress in thin film preparation has led to the pos- 
sibility of making high Q MgB2 microstrip structures op- 
erating at 20 to 30 K. This raises the question of how the 
nonlinearity of MgB2 compares with that of the cuprate 
superconductors 

Here, we study the intrinsic nonlinear microwave re- 
sponse of MgB2. By now it is well established that MgB2 
is a superconductor with two different superconducting 
gaps associated with different parts of the Fermi surface: 
a small gap (~ 2 meV) on the n band and a large gap 
(~ 7 meV) on the a band In a superconductor 

the intrinsic nonlinear response arises from the backflow 
of excited quasiparticles at finite temperatures. The to- 
tal current density j can be written as j = nev s — j qp 
where v s is the superfluid velocity and j qp the quasipar- 
ticle backflow. For a two band system like MgB2 the 
quasiparticle backflow consists of two contributions, one 
from each band with j qp — j qPtCr + j qP ,-*- In the clean 
limit these are given by 
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the Fermi velocity in each band, A Q (T) the gaps, / the 
Fermi function, N a (0) the partial densities of states, and 
(• • • } FS,a a Fermi surface average over band a. The term 
mvF,a ■ Vs describes the Doppler shift in energy of the 
quasiparticles with respect to the superflow. 

In MgB2 the Fermi surfaces of the two bands have 
significantly different topologies. In order to perform 
the Fermi surface averages in Q we use the band 
structure based model from 0, in which the a band 
is modeled as a distorted cylinder and the 7r band 
as a half torus. For the temperature dependence of 
the two gaps we are using the parameterization of 

Choi et. al A a (T) = A a (0)yjl- (T/T c ) 2 ' 9 and 

A W (T) = A 7r (0) y 1 - (T/Tc) 1 ' 8 with A CT (0) =6.8 meV 
and A^O) =1.8 meV. This was found to provide a rea- 
sonable fit to the anisotropic Eliashberg calculations of 

@ 

For currents flowing in the Boron plane, we expand 
Eq. up to third order in the current density. This 
gives the leading j 2 nonlinear change of the superfluid 
density HQ 



n s(T,j) = n s0 (T) 
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Here, the pair-breaking current densities of the two bands 
are given by j CiCe = eN a (0)vF, a A a (0). As discussed in 0, 
8] the nonlinear coefficients b a (T) determine the strength 
of the nonlinear response as a function of temperature 
and are given by 
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Here, a — n, a denotes the two bands in MgB2, vf.i 



Here, g a is a geometrical Fermi surface factor with 
5^=3.36 and g a =l. The linear part of the superfluid 
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FIG. 1: Temperature dependence of the total nonlinear co- 
efficient b(T) for MgB2 (solid line). The two contributions 
from the it band and the a band are shown as the dashed and 
dashed-dotted line, respectively. 
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FIG. 2: Temperature dependence of the intrinsic nonlinear 
coefficient b(T) for MgB2, for a d-wave superconductor with 
T c = 93 K (YBCO, dashed line), and a hypothetical BCS 
single gap superconductor with T c = 40 K (dashed-dotted 
line). 



density n s o(T) is given by 



n sQ {T) 



\ - n 



a— 7T,CT 



(4) 

where n a is the electron density of the a band. Using the 
known Fermi velocities and densities of states in MgB2, 
we find for the relative electron densities n„/n — 0.692 
and n a /n = 0.308 and the pair-breaking current densities 
j c>n = 3.32 • 10 8 A/cm 2 and j c , CT = 4.87 • 10 8 A/cm 2 . Inter- 
estingly, these values turn out to be close to the estimate 
of the pair-breaking current density of about 3T0 8 A/cm 2 
in YBCO, consistent with recent experiments 0. In 
Fig. n we show the temperature dependence of the total 

nonlinear coefficient b(T) — b 7T {T) J r-Sf :L b a {T) (solid line) 
along with the contributions from the it band (dashed) 
and the a band (dashed-dotted). Due to its smaller gap, 
the 7r band contribution dominates at low temperatures 
and leads to structure around 0.3T C . It is only at tem- 
peratures above about 0.75T C that the contribution from 
the cr band becomes significant. The total nonlinear co- 
efficient possesses a plateau-like region between 0.3 and 
0.7T C in which a reduction of the temperature does not 
improve the nonlinear response. 

In Fig.|5]we compare the nonlinear coefficient for MgB2 
(solid line) with the intrinsic response for a d-wave su- 
perconductor with T c — 93 K (YBCO, dashed line) 
and a hypothetical BCS single gap superconductor with 
T c = 40 K (dashed-dotted line) assuming for simplicity 
that the pair-breaking current densities are the same. For 
a ci-wave superconductor b(T) increases at low tempera- 
ture because of the gap nodes 0, Q • For this reason the 
nonlinear coefficient in MgB2 becomes smaller than the 
one for YBCO at about 27 K. However, due to the pres- 
ence of the small gap the nonlinear response in MgB2 is 



not as small as one would have expected for a single gap 
s-wave superconductor (dashed-dotted line). 

The foregoing analysis has been made in the clean limit 
without any impurity scattering. However, in most cur- 
rent MgB 2 films, the scattering rate as judged from the 
residual resistivity is larger than the two gaps [j| . For this 
reason, in the following we also want to discuss the dirty 
limit. In a two band superconductor there are in prin- 
ciple three different scattering rates: the two intraband 
scattering rates and and an interband scattering 
rate T^„. It has been shown that because of the different 
parity of the local orbitals making up the 7r and a bands, 
the interband scattering rate is much smaller than the 
intraband scattering rates 01- Thus, in the following we 
will neglect and keep only the intraband scattering 
rates. In the dirty limit the current densities j a in the 
two bands are then given by 
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where the renormalized Matsubara frequencies Cj n ,a and 
gaps A a are given by 
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and A a = T a 
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Expanding this expression up to third order in v s the 
sums over Matsubara frequencies can be done analyti- 
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FIG. 3: Temperature dependence of b(T) in the clean (solid 
line) and dirty limit. In the dirty limit b(T) depends on the 
relative scattering rates in the two bands of MgB2. Results 
are shown for Y^/Ya =0.7 (dotted), 1 (dashed-dotted), and 2 
(dashed). 
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FIG. 4: Dependence of b on ratio Y^/Ya at a fixed temper- 
ature of T = 0.5T C . For Y^ 3> r CT the nonlinear response is 
dominated by the a band, while for Y^ <C r CT the 7r band 
dominates. 



cally and we find 
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In the dirty limit the total nonlinear response coefficient 
b(T) = b„(T) + p^b a (T) only depends on the relative ra- 
tio of the two scattering rates Y^/Y^. In Fig. [21 we show 
its temperature dependence for T^/IV =0.7, 1, and 2 
along with the clean limit result. As can be seen, the non- 
linear response depends on the relative scattering rates 
in the two bands. In certain cases the coefficient b can 
even become smaller than in the clean limit. 

In order to have more insight into this unexpected 
behavior, in Fig. 0] we plot b at a fixed temperature, 
T = 0.5T C , as a function of the ratio Y^/Y^ on a double 
logarithmic scale. We find that b can vary by a factor 
of order 100 at this temperature. Qualitatively we can 
understand this behavior as follows: the current density 
is dominated by the band with the smaller scattering 
rate, because this band provides the highest conducting 
channel. If the it band scattering rate is smaller, the to- 
tal response is dominated by the small gap leading to a 
larger nonlinear response at finite temperature. If, how- 
ever, the a band scattering rate is smaller, the nonlinear 
response is dominated by the large gap giving a smaller 
nonlinear response. 

This analysis tells us that optimization of material 
properties in MgB2 should aim at a higher scattering 



rate in the n band in order to suppress the contribution 
of the small gap relative to that of the large gap in the 
a band. This could be achieved by substitutional doping 
at the Mg site, for example with Aluminum 10]. Our 
analysis also shows that clean MgB2 does not necessarily 
provide the lowest nonlinear response. 

As has been shown in |2j the intermodulation power 
emitted by a microstrip resonator also depends on the 
penetration depth of the material. A shorter penetra- 
tion depth leads to an increased current density at the 
edges of the resonator, which increases the intermodula- 
tion power. According to Refs. and [j| the intermod- 
ulation power Pimd scales like 



Pimd oc (A£) 2 b 2 (T) 



(9) 



where AC is a nonlinear coefficient depending on both 
the geometry of the resonator and the penetration depth. 
According to band structure calculations, the clean limit 
zero temperature (London) penetration depth in MgB2 
is expected to be near Ai(0) =40 nm [[J. Actual values 
vary between 60 and 200 nm depending on film qual- 
ity Comparing MgB 2 with the intrinsic d-wave 
response of YBCO at a temperature of 20 K, we take 
A(T = 20K) « 100 nm for MgB 2 and A(T = 20^) w 
160 nm for YBCO. Assuming a typical film thickness of 
t = 400 nm we find for the microstrip geometry con- 
sidered in Ref. 2] AC(t/X = 4) = 0.173 (MgB 2 ) and 
AC(t/X = 2.5) = 0.124 (YBCO). Taking b(T = 20K) = 
0.563 for MgB 2 and b(T = 20K) = 0.855 for YBCO from 
Fig. [21 and using Eq. © this means that the intermodula- 
tion power in YBCO at this temperature would be larger 
by only a factor of 1.2, i.e. a comparable nonlinear re- 
sponse for both materials. However, due to their ceramic 
nature, weak links play a much larger role in the high-T c 
cuprates than in MgB 2 . For this reason we expect that 
the intrinsic nonlinear response will be much easier to 
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achieve in MgB2 films compared with cuprate films. 
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